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The rapid progress of heavy ion collision experiments 
urge us towards quantitatively reliable determination of 
these properties. Lattice QCD is a powerful basis to 
investigate thermodynamic properties of the quark matter 
around the deconfinement transition temperature directly 
from the first principles of QCD. There were quite a few 
big advances in finite temperature QCD on the lattice 
made in the last years. 

Because the s quark mass is comparable with the QCD 
scale Aqcd and thus the transition temperature, it is im- 
portant to incorporate the s quark dynamically for quan- 
titative predictions. Large-scale simulations in 2+1 fla- 
vor QCD with various improved staggered quark ac- 
tions have now been started to produce results for var- 
ious thermodynamic quantities which are extrapolated to 
the continuum limit at around physical quark masses. At 
the same time, the theoretical uneasiness with staggered 
quarks motivated several groups to accelerate studies 
with Wilson-type quarks and with lattice chiral quarks. 
I review these recent developments in finite temperature 
lattice QCD, concentrating on the topics of the phase 
structure and the equation of state. See [1] for more dis- 
cussions. 



TRANSITION TEMPERATURE 

The transition temperature T c is one of the most impor- 
tant quantities for experimental investigations of QGP. 
Estimation of T c in 2+1 flavor QCD has been made based 
on large-scale simulations using various improved stag- 
gered quarks. However, there has been a sizable discrep- 
ancy in T c among groups for more than five years. 

The discrepancy. In 2005, the MILC Collaboration 
obtained T c = 169(12) (4) MeV in the chiral and con- 
tinuum limits from a study of the chiral susceptibility 
using asqtad quarks [2] on N t — 4-8 lattices [3], where 
the scale was set by r\ and the 0(4) critical exponents 



were adopted in the chiral extrapolation. In 2006, based 
on a simulation using stout quarks [4] on N t = 6-10 lat- 
tices, the Wuppertal-Budapest Collaboration published 
T c = 15 1 (3) (3) MeV from the susceptibility of the chiral 
condensate and 175 (2) (4) MeV from the s quark num- 
ber susceptibility [5]. These are the results extrapolated 
to the physical point and to the continuum limit, and 
the scale was set by fg. In the same year, the BNL- 
Bielefeld Collaboration published T c = 192(7) (4) MeV 
at the physical point in the continuum limit, based on a 
study using p4 quarks [6] on N t — 4 and 6 lattices [7]. 
This T c is an average of the values from the Polyakov- 
loop susceptibility and the chiral susceptibility, and the 
difference between them is included in the systematic 
error. Later, the HotQCD Collaboration (a merger of 
the BNL-Bielefeld and MILC Collaborations) compared 
p4 and asqtad actions and found that the two staggered 
quark actions lead to roughly consistent T c on finite lat- 
tices, while it shifts towards smaller values when the lat- 
tice spacing is decreased and when the light quark mass 
is decreased [8, 9, 10]. 

Because the transition is considered to be an analytic 
crossover around the physical point (see e.g. [11]), the 
value of T c may depend on the choice of observable to 
define it. T c should depend on the details of the analyses, 
such as the scale convention, too. However, it was found 
that the discrepancy exists even when we adopt the same 
observable and the same scale convention. 

Resolution. The Wuppertal-Budapest Collaboration 
extended the study to finer lattices (N t up to 16) and 
concluded that their values remain essentially unchanged 
[12, 13]. Their final values are T c = 147(2)(3) MeV 
from the susceptibility of the chiral condensate and T c = 
165(2) (4) MeV from the s quark number susceptibil- 
ity [13]. At the Lattice 2010 conference, the HotQCD 
Collaboration presented a preliminary result T c = 164(6) 
MeV at the physical point in the continuum limit, based 
on a study of the disconnected chiral susceptibility using 




FIGURE 1. Chiral observables in 2+1 flavor QCD with improved staggered quarks. (Left) Subtracted chiral condensate 
(difference of scaled light quark and s quark condensates to remove a divergent renormalization factor) from the HISQ, p4, asqtad 
and stout actions [14]. (Center) Comparison of HISQ and asqtad results for the disconnected part of the chiral susceptibility [14]. 
In these panels, the HISQ, p4 and asqtad data are obtained at the bare quark mass ratio m llc i /m s = 0.05 with m s around the physical 
value, using a scale fixed by r$ . The stout data, which is obtained at the physical point using a fx scale [13], is shifted in T to 
correct the difference in the scale setting. Note that this procedure causes a slight deviation from the physical point at finite lattice 
spacings. (Right) Subtracted chiral condensate with stout and asqtad quarks at N, = 8 [13]. The open rectangles are results of stout 
quark action at mJ NG = 414 MeV [13], and the red curve is the result of asqtad quark action at m^ NG = 220 MeV [9]. The stout 
quark mass for open rectangle is adjusted to reproduce m^ MS = 587 MeV of the asqtad data at T ~ 135 MeV on an N, = 8 lattice. 



asqtad quarks on lattices up to N t = 12 [14]. In parallel, 
they presented first results of their new project [14, 15] 
adopting the "highly improved staggered quark (HISQ)" 
action [16]. Comparison of these actions is shown in left 
and center panels of Fig. 1. Consulting new results, we 
find that the discrepancy in T c is now mostly removed: 
All the actions converge towards T c w 145-165 MeV. 

Staggered quarks and the taste violation. What was 
the origin of the discrepancy? At this point, I have to 
comment on several caveats with staggered-type lattice 
quarks. A naive lattice discretization of the Dirac ac- 
tion leads to the species doubling problem. To avoid 
this, various lattice fermions have been proposed, in- 
cluding the Wilson-type, the staggered-type and lattice 
chiral fermions. Among them, most quantitative out- 
puts of finite temperature QCD have been produced with 
staggered-type lattice quarks because of their relatively 
cheap simulation costs. 

On the other hand, the staggered-type quarks lead to 
four copies ("tastes") of fermions for each flavor in the 
continuum limit. To remove unwanted three tastes, we 
usually adopt the "fourth root procedure", i.e. we re- 
place the quark determinant by its fourth root in the path- 
integration procedures. Since this makes the theory non- 
local and non-unitary, many theoretical issues emerge 
[17], as intensively discussed in the Confinement 8 work- 
shop. In particular, the universality arguments become 
fragile due to non-local interactions, and thus, e.g. we are 
not fully sure if the rooted staggered quarks have the cor- 
rect continuum limit, though, so far, no apparent obsta- 
cles about the validity of the rooting procedure seem to 
exist, provided that the continuum extrapolation is done 
prior to the chiral extrapolation [17, 18]. 



Even when we assume that the staggered quarks do 
have the correct continuum limit of QCD with desired 
number of flavors, a couple of issues remain on finite 
lattices. One of them is the problem of taste violation. 
The symmetry among tastes in the continuum limit is ex- 
plicitly broken at finite lattice spacings. This introduces 
a systematic uncertainty in the identification of flavors, 
and we observe different masses for a hadron depending 
on the combination of tastes. For the pion, one combina- 
tion of tastes, the pseudo Nambu-Goldstone (pNG) pion, 
leads to the vanishing mass in the chiral limit, but other 
15 combinations remain heavy. Although the pNG pion 
is customary treated as the pion, heavier pions do con- 
tribute to loop diagrams and thus induce lattice artifacts. 
To reduce the lattice artifacts, various improved actions, 
including the asqtad, p4, stout and HISQ actions, have 
been proposed. 

Because loop diagrams are the leading contributions 
in many thermodynamic quantities, a good control of the 
taste violation is important in finite temperature QCD. 
In the studies of T c , the main difference among the 
staggered-type quark actions seems to be the magnitude 
of the remaining taste violation: When we adjust the 
quark mass parameter to have w^ NG « 135 MeV at lat- 
tice spacings corresponding to T ss 170 MeV on TV, = 8 
lattices, the masses for the heavy pions are around 400- 
600 (asqtad), 300-500 (stout) and 200-400 MeV (HISQ), 
respectively. The p4 is slightly worse than the asqtad. 
Therefore, the p4 and asqtad suffers from severer lattice 
artifacts than the stout and HISQ. 

The overall effects of heavy pions in the loop diagrams 
may be measured by a root mean squared (RMS) mass of 
pions m^ MS . In the right panel of Fig. 1, the asqtad data 
is shown to be well reproduced by the stout quark action 
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FIGURE 2. Order of the finite temperature transition in 2+1 
flavor QCD as a function of the degenerate u and d quark mass 
m uc i and the s quark mass m s . (Left) The standard scenario with 
the second order chiral transition for two-flavor QCD. (Right) 
An alternative scenario when the two-flavor chiral transition is 
first order. 



by adjusting m^ MS [13]. This explains the larger T c with 
the p4 and asqtad actions at A', < 8. 

While the major part of the discrepancy is explained 
by the taste violation, the sensitivity of the determination 
of T c on the operator, observed with the stout quark, is 
not confirmed by other quarks yet. 

In all of these studies, the physical point is identified 
by the pNG pion mass. The success of the RMS pion 
mass suggests that it is more appropriate to consult the 
RMS hadron masses to judge the location of the physical 
point for thermodynamic quantities. This will make the 
physical point closer to the chiral limit and the values of 
T c even smaller on finite lattices. 

Other lattice quarks. It is highly desirable to confirm 
the results with other types of lattice quarks whose the- 
oretical basis is solid. With Wilson-type lattice quarks, 
the correct continuum limit is guaranteed for any number 
of flavors, at the price of explicit violation of the chiral 
symmetry at finite lattice spacings and more computa- 
tional costs. Due to the higher computational costs, the 
status of the studies with Wilson-type quarks is quanti- 
tatively much behind at present. Recent updates of T c in 
two-flavor QCD are given in Refs. [19, 20, 21]. 

The most attractive but the most expensive way is to 
adopt a chiral lattice quark action, such as the domain- 
wall and overlap. The first estimate of T c with 2+1 flavors 
of domain-wall quarks is made by the HotQCD Collab- 
oration [22]. The JLQCD collaboration has started finite 
temperature studies of QCD using an overlap quark ac- 
tion at a fixed topology [23]. 



PHASE STRUCTURE 

Figure 2 summarizes the current wisdom about the order 
of the finite temperature transition in 2+1 flavor QCD as 



a function of the quark masses, based on the studies on 
the lattice and with effective models. When the u, d and 
s quarks are all sufficiently light or sufficiently heavy, 
we expect the transition to be of first order. At interme- 
diate values of quark masses, the "transition" will be an 
analytic crossover. On the boundaries of the first order 
regions, we expect second order transitions. Associated 
critical scaling around there will have characteristics uni- 
versal to the corresponding effective models. 

The nature of the transition in the chiral limit of two- 
flavor QCD (the upper left edge of the figure) has signif- 
icant implications for the nature of the transition at the 
physical point too. The transition is predicted to be either 
second order in the universality class of the 0(4) Heisen- 
berg model, or first order if the anomaly is weak around 
the transition point [24] . 

The left panel of Fig. 2 summarizes the standard sce- 
nario in which the chiral transition of two-flavor QCD is 
second order. In this case, because the chiral transition 
of three-flavor QCD is of first order, we have a tricritical 
point on the left edge of the figure (m u[ i = 0) where the 
order of the transition changes from the second order to 
the first order. Depending on the location of the tricrit- 
ical point relative to the physical point, the universality 
class dominating the parameter dependence around the 
physical point will be different. 

The right panel of Fig. 2 shows an alternative scenario 
in which the chiral transition of two-flavor QCD is first 
order. In this case, we have no tricritical point and thus no 
regions for the 0(4) universality class. A distinction be- 
tween the two scenarios is important for studies at finite 
densities too. Although the majority view the standard 
scenario as more probable, the nature of the two-flavor 
chiral transition was not fully fixed. Here, scaling is a 
powerful tool to discriminate the nature of the transition. 

0(4) scaling with Wilson-type quarks. With Wilson- 
type quarks, the 0(4) scaling has been observed for var- 
ious actions since the early years of the investigations. 
Although Wilson-type quarks explicitly violate the chi- 
ral symmetry at finite lattice spacings, when the theory is 
correctly renormalized such that the chiral violations are 
subtracted out and the chiral Ward identities are satisfied, 
the resulting subtracted observables should show the cor- 
rect chiral behaviors near the continuum limit [25]. 

Adopting a chiral condensate defined by the axial vec- 
tor Ward identity [26], the QCD data are shown to sat- 
isfy the 0(4) scaling, with both unimproved and clover- 
improved Wilson quark actions [26, 27, 28]. See the left 
panel of Fig. 3. A recent study on finer lattices with 
lighter quark masses supports the 0(4) behavior too [20]. 
These results, together with inconsistencies with first or- 
der scaling ansatze, suggest that the transition is contin- 
uous in the chiral limit. However, the quark masses stud- 
ied are heavy yet (»% > 400 MeV). Confirmation with 
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FIGURE 3. (Left) 0(4) scaling test for the subtracted chi- 
ral condensate in two-flavor QCD using the clover-improved 
Wilson quark action and the RG-improved Iwasaki gauge ac- 
tion on an N t = 4 lattice by the CP-PACS Collaboration[27]. 
The dashed curve is the scaling function obtained in an 0(4) 
spin model. Similar result supporting the 0(4) scaling was 
obtained with the unimproved Wilson quark action and the 
Iwasaki gauge action too [26]. (Right) Rescaled chiral con- 
densate Mq = m s { iffy) U d/T A in 2+1 flavor QCD using the 
p4-improved staggered quark action and tree-level Symanzik 
gauge action on an N t = 4 lattice by the BNL-Bielefeld Col- 
laboration [34]. The QCD data are fitted to the 0(2) scaling 
function with correction terms modeling a deviation from the 
scaling. 



lighter quarks on a fine lattice is indispensable. Studies 
towards this direction are in progress (see e.g. [29]). 

Chiral scaling with staggered-type quarks. Scal- 
ing studies with staggered-type quarks have additional 
caveats. With the fourth root procedure, the symmetry 
of the system at finite lattice spacings is identical to that 
of the original four-taste theory for any Np. In the chi- 
ral limit, we thus have the 0(2) symmetry. This means 
that, when the chiral transition is of second order, we 
should expect scaling properties in the universality class 
of 0(2).' The 0(4) scaling may appear only when we 
take the continuum limit prior to the chiral extrapolation. 
We may, however, expect that when the chiral transition 
is continuous on finite lattices, it will remain so in the 
continuum limit. In this sense, examination of 0(2) is im- 
portant, though the 0(2) scaling is numerically not easy 
to discriminate from the 0(4) scaling. 

The results for the chiral transition with two flavors 
of staggered quarks have been quite confusing. Although 
the early studies suggest a continuous chiral transition, 
the scaling observed was neither 0(2) nor 0(4) [30, 31, 
32]. In contradiction to these studies, the Genova-Pisa 
group reported indications of a first order transition on 
anN t =4 lattice [33]. 

All of these studies are made using the unimproved 



1 In the studies of the 0(2) scaling with staggered-type quarks, it is 
assumed that the non-locality of rooted staggered quark actions does 
not invalidate the universality arguments at finite lattice spacings too, 
though no theoretical argument supporting this assumption is known. 



staggered quark action and the unimproved plaquette 
gauge action, however. As discussed in the previous sec- 
tion, large taste-violation with these actions will make 
the quark masses effectively much heavier than those 
naively estimated by the pNG pion mass. 

This year, a new study with an improved staggered 
quark action was presented: Using p4 quarks, the BNL- 
Bielefeld Collaboration performed a scaling test in 2+1 
flavor QCD on an N t — 4 lattice for the quark mass range 
corresponding to m^ NG w 150 MeV down to 75 MeV 
[34], which is much lighter than the previous studies. 
For the first time, the 0(2) scaling for staggered-type 
quarks has been observed, as shown in the right panel of 
Fig. 3. At the same time, they found substantial deviation 
from the 0(2) scaling for m uc i/m s > 0.05. Together with 
the adoption of improved actions, this deviation explains 
the failures in previous studies with much heavier quark 
masses. A preliminary result obtained on a finer N t = 8 
lattice was shown to be similar [35]. 

A big impact of this study is the 0(2) dominance 
around the physical point. First, this suggests that the 
chiral transition in two-flavor QCD is second order. Sec- 
ond, the s quark mass for the tricritical point may be 
smaller than the physical value. In Fig. 2, I have there- 
fore put the physical point slightly higher than the pre- 
vious plots. When the parameter-dependence around the 
physical point is confirmed to be dominated by the 0(4) 
scaling in the continuum limit, the studies of QGP will 
be largely accelerated. 

Near the heavy quark limit. The upper right corner 
of Fig. 2 corresponds to the heavy quark limit where the 
deconfinement transition of SU(3) pure gauge theory is 
described by the first order transition of an effective Z(3) 
Potts model [36]. Dynamical heavy quarks affect as an 
external magnetic field to the effective spins and break 
the Z(3) symmetry explicitly. 

The WHOT-QCD Collaboration studied the fate of 
the deconfinement transition when the quark mass is 
decreased from infinity [37]. Adopting a method using 
probability distribution of observables combined with 
a reweighting technique, they succeeded in calculating 
the effective potential, defined as the logarithm of the 
probability distribution for the plaquette, for a wide range 
of the plaquette expectation value. Close to the heavy 
quark mass limit, the slope of the effective potential 
shows the shape like the letter "S" that is typical for a 
first order transition. The "S"-shape is shown to become 
weaker and eventually invisible when m q is decreased. 
Determining the location of the critical point terminating 
the first order transition, they obtained the phase diagram 
for 2+1 flavor QCD. As naively expected, the endpoint 
monotonically shifts to lighter quark mass as we increase 
the number of flavors. 



EQUATION OF STATE 

Calculation of the EOS in full QCD is demanding. Be- 
sides the basic information such as the lattice scale, the 
non-perturbative beta functions and the line of constant 
physics (LCP), we need the expectation values at T = 
to renormalize finite temperature observables at each 
simulation point. Conventionally, we vary the coupling 
parameters along a LCP with a fixed temporal lattice size 
Nt to vary the temperature (fixed N t approach). To calcu- 
late the EOS non-perturbatively, the integration method 
is adopted [38]. These require a systematic scan in a wide 
range of parameters on both zero and finite temperature 
lattices with sufficiently large spatial sizes. For a contin- 
uum extrapolation, we repeat the calculation at several 
values ofN,. 

EOS with 2+1 flavors of improved staggered quarks. 
The BNL-Bielefeld Collaboration published the EOS 
"with physical quark masses", based on their simulations 
using p4 quarks on an N t = 8 lattice [10]. Using a scale 
set by ro, their simulation point at the bare quark mass 
ratio m ul j/m s = 0.05 corresponds to w^ NG = 154 MeV. 
The Wuppertal-Budapest Collaboration calculated the 
EOS with stout quarks at the physical point on N t = 6- 
10 lattices (and a few points at N t = 12) [39]. They 
also attempted to calculate the EOS up to T = 1000 
MeV using LCP determined by a step-scaling method. 
At Lattice 2010, the HotQCD Collaboration presented 
preliminary results of the EOS with HISQ and asqtad 
quarks at m uc i/m s = 0.05 on N t = 8 lattices [14]. 
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FIGURE 4. Trace anomaly in 2+1 flavor QCD. The scale is 
set by ro, except for the stout data for which the scale is set 
by fx. (Left) Comparison of HSQ, asqtad and p4 quarks at 
m ud/ m s = 0.05 on N, = 8 lattices [14]. (Right) Comparison of 
stout, asqtad and p4 quarks on lattices V ; > 8 [39]. The stout 
data is obtained at the physical point (m uc j/m s ~ 0.035) while 
asqtad and p4 data are for m uc j/m s = 0.05. The stout values are 
corrected by a tree-level improvement factor. 

The results for the trace anomaly from these studies 
are plotted in Fig. 4. In the left panel, a comparison is 
made at m uc i/m s = 0.05 for HSQ, asqtad and p4 quarks. 
In the right panel, the stout data is compared with p4 
and asqtad results. While the location of the transition is 
approximately consistent with each other, the magnitude 
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FIGURE 5. EOS in 2+1 flavor QCD. (Left) p4 results at 
mP NG & 154 MeV {m ud /m s = 0.05) obtained on an N, = 8 
lattice using the ro scale [10]. Results at m^ NG = 220 MeV 
(m uc i/m s = 0.1) are compared. (Right) Pressure with stout 
quarks using the fg scale [39]. 



of the trace anomaly shows clear dependence on the 
action adopted. In particular, we note that the peak of 
the trace anomaly around T ~ 200 MeV becomes lower 
with the adoption of a more improved action. 

Using the integration method, the EOS around the 
physical point is obtained, as shown in Fig. 5. Note 
that, in these studies, the physical point was identified 
by the pNG pion mass, and the chiral and continuum 
extrapolations are not done yet. 

EOS with 2+1 flavors of improved Wilson quarks. 
About ten years ago, the CP-PACS Collaboration per- 
formed a systematic calculation of the EOS in two- 
flavor QCD with clover-improved Wilson quarks cou- 
pled to RG-improed Iwasaki glue [27, 40]. Extension 
of the study to 2+1 flavor QCD was not straightforward 
due to the heavy computational demand for Wilson-type 
quarks. 

To overcome the problem, the WHOT-QCD Collabo- 
ration developed a new approach: a fixed scale approach 
[41] in which T is varied through a variation of N t with 
the coupling parameters fixed. The conventional integra- 
tion method [38] is inapplicable in the fixed-scale ap- 
proach. Therefore, they developed an alternative inte- 
gration method, the T-integration method [41]. With the 
fixed scale approach, the cost for zero-temperature simu- 
lations can be largely reduced, because a common zero- 
temperature result can be used to the renormalization 
at all temperatures. We may even borrow high statistics 
zero-temperature configurations publicly available on the 
International Lattice Data Grid (ILDG). The fixed-scale 
approach is complementary to the fixed N t approach. To- 
wards the high T limit, the fixed-scale approach suffers 
from lattice artifacts due to small N t , while the fixed N t 
approach is applicable for quantities insensitive to the 
system volume. In the low T region, the fixed N t ap- 
proach suffers from coarseness of the lattice, while the 
fixed-scale approach keeps the lattice spacing. 

Adopting the fixed-scale approach, the WHOT-QCD 
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FIGURE 6. EOS in 2+1 flavor QCD with clover-improved 
Wilson quarks adopting the fixed-scale approach and the T- 
integration method [42]. 



Collaboration attempted the first calculation of the EOS 
in 2+1 flavor QCD with Wilson-type quarks [42]. They 
borrowed zero-temperature configurations using the non- 
perturbatively improved clover quark action [43], and 
performed finite temperature simulations on N t = 4-16 
lattices, which corresponds to T s» 170-700 MeV. Their 
preliminary result for the EOS is shown in Fig. 6. The 
peak height of the trace anomaly, obtained on N t ~ 
14 lattices, is roughly consistent with those shown in 
Fig. 4 with highly improved staggered quarks, though 
the statistical error is large. Further statistics are needed 
to reduce the errors. 



SUMMARY 

There have been important advances in finite tempera- 
ture QCD on the lattice in the past years. Large-scale 
simulations of 2+1 flavor QCD with various improved 
staggered quark actions clarified the origin of the recent 
discrepancy in the transition temperature. From this, we 
have learned the importance of the effect of taste viola- 
tion in the studies with staggered-type quarks. The 0(2) 
chiral scaling for staggered-type quarks was observed 
by reducing the light quark mass down to a quite low 
value. At the same time, the theoretical uneasiness with 
staggered-type lattice quarks motivated several groups to 
accelerate studies adopting Wilson-type and lattice chi- 
ral quark actions. Calculation of the EOS in 2+1 flavor 
QCD has started with improved Wilson quarks. Further- 
more, simulations with domain-wall quarks have entered 
a level of quantitative investigations. 
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